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ABSTRACT 

We present low-resolution (R~ 850) spectra for 67 asymptotic giant branch, horizontal branch and red 
giant branch (RGB) stars in the low-metallicity globular cluster NGC 5466, taken with the VIRUS-P 
integral-field spectrograph at the 2.7-m Harlan J. Smith telescope at McDonald Observatory. Sixty-six 
stars are confirmed, and one rejected, as cluster members based on radial velocity, which we measure 
to an accuracy of 16 kms _1 via template-matching techniques. CN and CH band strengths have been 
measured for 29 RGB and AGB stars in NGC 5466, and the band strength indices measured from 
VIRUS-P data show close agreement with those measured from Keck/LRIS spectra previously taken 
of five of our target stars. We also determine carbon abundances from comparisons with synthetic 
spectra. The RGB stars in our data set cover a range in absolute V magnitude from +2 to —3, 
which permits us to study the rate of carbon depletion on the giant branch as well as the point of its 
onset. The data show a clear decline in carbon abundance with rising luminosity above the luminosity 
function "bump" on the giant branch, and also a subdued range in CN band strength, suggesting 
ongoing internal mixing in individual stars but minor or no primordial star-to-star variation in light- 
element abundances. 

Subject headings: globular clusters: general - globular clusters: individual (NGC 5466) - stars: abun- 
dances 



1. INTRODUCTION 

In all Galactic globular clusters with moderate metal- 
licity ([Fc/H] ~ —1-5), there are star-to-star variations 
in the strength of the 3883A CN absorption band at any 
given evolutionary phase. These CN band strengths an- 
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ticorrelate with CH band strengths, so that stars with 
weak CH bands (and therefore relatively low carbon 
abundances) have strong CN bands (and therefore rel- 
atively high nitrogen abundances) . Band-strength varia- 
tions were first observed in spectra of red giants in globu- 
lar clusters (e.g., Norris & Freeman 1979, Suntzeff 1981), 
but are also found in subgiants and on the main sequence 
(e.g., Briley et al. 2004), indicating that they are a result 
of primordial enrichment rather than an evolutionary ef- 
fect. 

Variations in the broad CN and CH features are just 
the most readily observed part of a larger light-element 
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abundance pattern, a division of globular cluster stars 
into one group with typical Population II abundances 
and another that is relatively enhanced in nitrogen, 
sodium, and aluminum, and depleted in carbon, oxy- 
gen, and magnesium. These abundance divisions are 
usually studied in the anticorrelated pairs of carbon and 
nitrogen or oxygen and sodium (see, e.g., Kraft 1994; 
Gratton et al. 2004; Carretta et al. 2009). The 
magnesium-aluminum anticorrelation is more murky: al- 
though some researchers find Mg-Al anticorrelations in 
individual globular clusters (e.g., Shetrone 1996; Grat- 
ton et al. 2001), the abundance distributions are not as 
clearly bimodal as in the case of nitrogen, and the ori- 
gin of Mg-poor, Al-rich stars is more difficult to fit into 
the general framework of enrichment by moderate-mass 
AGB stars (e.g., D'Ercole et al. 2008). 

This paper presents a study of the behavior of CN and 
CH band strengths and [C/Fe] abundances in red giant 
stars in the low-metallicity globular cluster NGC 5466, 
which has [Fe/H] ~ 2.2 (Harris 1996, 2003 revision). 
Studies based on CN and CH band strength are more dif- 
ficult in low-metallicity clusters, where the CN variation 
visible in the spectra can be quite small despite variations 
in [N/Fe] abundance as large as in higher-metallicity 
globular clusters (e.g., M53, Martell et al. 2008a; M55, 
Briley et al. 1993). The globular clusters with metallic- 
ities similar to NGC 5466 in which light-element abun- 
dance inhomogeneities have been most extensively stud- 
ied are M92 and M15. The red giants in both objects 
exhibit a progressive decline in [C/Fe] with advancing 
evolution above the magnitude level of the horizontal 
branch (HB; Carbon et al. 1982; Trefzger et al. 1983; 
Langer et al. 1986; Bellman et al. 2001). Giants in both 
clusters typically exhibit weak A3883 CN bands (Carbon 
et al. 1982; Trefzger et al. 1983), and do not exhibit 
bimodal CN distributions like those in more metal-rich 
clusters although a handful of stars with enhanced CN 
have been discovered in M15 (Langer, Suntzeff, & Kraft 
1992; Lee 2000). Despite their lack of CN-strong giants, 
both M92 and M15 exhibit anticorrelated O-N, O-Na or 
O-Al variations, or correlated N-Na, of the type that are 
commonplace in globular clusters (Norris & Pilachowski 
1985; Sneden et al. 1991, 1997; Shetrone 1996). In ad- 
dition, Cohen, Briley & Stetson (2005) discovered that 
star-to-star differences in [C/Fe] and [N/Fe] abundances 
exist among stars near the base of the red giant branch 
(RGB) in M15, and the C and N abundances tend to be 
anticorrelated. These efforts have shown that for very 
metal-poor stars the use of CN to determine which stars 
carry the anticorrelated light-element abundance pattern 
can be problematic, and that anticorrelated variations 
in light-element abundances do exist in low-metallicity 
globular clusters despite their lack of obvious CN band 
strength variation. 

If NGC 5466 does contain the same primordial varia- 
tions in light-element abundances as are found in higher- 
metallicity globular clusters, the resulting range in CN 
and CH band strength will be muted by the low overall 
metallicity of the cluster, perhaps to a degree where it is 
difficult to distinguish a bimodal but closely spaced CN 
band strength distribution from a broad but unimodal 
distribution. If NGC 5466 does not contain primordial 
abundance variations, that would imply that primordial 
enrichment is not a universal process in low-metallicity 



globular clusters the way it is in higher-metallicity clus- 
ters. A metallicity limit on primordial light-element en- 
richment offers some insight into the primordial enrich- 
ment process, and into the larger and still-developing 
picture of chemical complexity within individual globu- 
lar clusters (e.g., Piotto 2008). 

2. OBSERVATIONS AND REDUCTION 

We determined photometry for members of NGC 5466 
based on CCD frames obtained by HEB with the Kitt 
Peak National Observatory 0.9 m telescope on 1997 May 
9. The T2KA chip at the Ritchey-Chretien focus pro- 
vides a 23' x 23' field, which encloses most of the clus- 
ter, which has a tidal radius of 34', according to Harris 
(1996). We used uBVI filters, as defined by Bond (2005), 
with exposure times of 600, 45, 45, and 60 s, respectively, 
under photometric conditions. We calibrated the pho- 
tometry to the network of uBVI standard stars estab- 
lished by Siegel & Bond (2005). Data were reduced using 
the IRAF CCDPROC pipeline and photometry was mea- 
sured with DAOPHOT / ALLSTAR (Stetson 1987, 1994). 
DAOGROW (Stetson 1990) was used to perform curve- 
of-growth fitting for aperture correction on both program 
and standard stars. The raw photometry was calibrated 
using the iterative matrix inversion technique described 
in Siegel et al. (2002) to translate the photometry to the 
standard system of Landolt (1992) and Siegel & Bond 
(2005). Table Q] lists identification numbers, positions, 
and photometry for the stars observed with VIRUS-P. 
The star identification numbers and coordinates are from 
an unpublished photometric catalog developed by M.H.S. 
from the analysis described above. The typical errors on 
the B, V, and I photometry are 0.030, 0.022 and 0.032 
mag, respectively. 

On 2009 March 3-5 and July 28-29, we observed four 
fields in NGC 5466 for a total of 22 hr with the visi- 
ble integral-field replicable unit spectrograph prototype 
(VIRUS-P, Hill et al. 2008) on the 2.7m Harlan J. Smith 
Telescope at McDonald Observatory. VIRUS-P is a spec- 
trograph with an integral field unit (IFU) of 246 fibers 
in a hexagonal, one third fill factor, close pack pattern. 
The fibers have a projected diameter of 4.1" and have 
centers astrometrically calibrated to 0.7" relative to a 
fixed, offset guiding camera via open cluster observa- 
tions. The field of view for the IFU is 2.89 square arcmin 
and with a set of three 7.14" dithers the filling factor is 
100%. The VP1 grating was used, giving an instrumen- 
tal FWHM between 5 and 6A over the wavelength range 
of 3500 — 5800A. Arc lamp calibration produced a wave- 
length solution with an rms accuracy of 0.05A. Our pro- 
cedure during observing was to dither the telescope and 
take three 1200 s exposures per field to fill in the spaces 
between the IFU fibers. Three of the fields had four sets 
of dithers, but the July observing run was less successful 
due to weather and only a single dither set was taken for 
the final field. 

A custom software pipeline developed for the instru- 
ment (Adams 2010, in preparation) was used in all re- 
ductions. The notable features are that no interpolation 
was performed on the data to avoid correlated noise and 
that background subtraction was done between differ- 
ently sampled fibers by fitting B-spline models (Dierckx 
1993) in a running boxcar of 31 chip-adjacent fibers in a 
manner similar to optimal longslit sky subtraction meth- 
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Fig. 1. — Color magnitude diagram for the photometric sample 
of NGC 5466. Filled squares represent the spectroscopic sample, 
and open squares show stars with high enough S/N for spectro- 
scopic analysis but the possibility of spectroscopic contamination 
from a nearby star. The open triangles denote HB stars and red 
giants with insufficient S/N for spectral analysis. The single upside 
down triangle represents star 1980, a strong-lined RV non-member. 
Small points represent stars not observed as part of this VIRUS- 
P campaign. Filled squares outlined by larger symbols represent 
AGB stars (outlined with circles) and a probable CH star (outlined 
with a star) . 

ods (Kelson 2003). The NGC 5466 pointing was sparse 
enough to allow self-sky-subtraction. The spectra were 
flux calibrated by observing a single flux standard on 
each night using the same dithering pattern used for the 
science targets. The flux solution was then applied to 
every fiber in every exposure for that night. Based on 
many tens of flux standards taken during other VIRUS- 
P projects the absolute flux calibration precision is 8.5% 
and wavelength independent for full dither sets. The 
data reduction code does not correct for differential at- 
mospheric refraction (DAR) because there is very little 
impact from DAR with these large fibers when they are 
summed over full dither sets. Fiber spectra were then 
combined using the IRAF task scombine to form a single 
spectrum for each star. For some stars, this was done 
by combining a single fiber from all visits (usually four) 
to a specific dither position; for other stars, those falling 
near the edge of a fiber, it may have required combin- 
ing different fibers on different dithers (perhaps 8 or f 2 
spectra). 

Figure [1] is a color-magnitude diagram for NGC 5466, 
with different symbols denoting various parts of our 
data set: stars not observed in this particular VIRUS-P 
pointing are small points, and the spectroscopic sample 
(unblended asymptotic giant branch (AGB) and RCB 
stars with high enough signal-to-noise for index mea- 
surements) are shown as filled squares. Our standard 
for a spectroscopically unblended star is that there are 
no other stars within 5 arcsec down to 3 mag below 
the sample star's B magnitude. Open squares repre- 
sent stars with sufficient signal-to-noise (S/N) for spec- 
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Fig. 2. — Four spectra from our spectroscopic sample as examples 
of the characteristics of VIRUS-P spectra. The top three spectra 
were chosen to span the range of absolute magnitude in our data 
set. The lowest spectrum, that of star 1839, has very strong CN and 
CH bands. 1839 also sits to the red of the RGB in the (V, B-V) 
CMD, and is probably a CH star. 

troscopic analysis but possible spectroscopic blends, and 
open triangles represent HB stars and RGB stars with 
low S/N. One strong-lined radial velocity nonmember 
(star 1980) is shown as an inverted triangle. The two 
filled squares outlined with large circles are stars 1398 
and 2483, and are likely to be AGB stars based on their 
positions to the blue of the giant branch. The filled 
square outlined with a large star is star 1839, which has 
strong CN and CH features and is likely a CH star. Al- 
though 1839 lies on the giant branch in the (V, V — I) 
color-magnitude diagram, the strong absorption in the 
UV, due to extra CN opacity and the Bond-Neff effect 
(Bond & Neff 1969), makes its (B - V) color distinctly 
red relative to the giant branch. The distance modulus 
and reddening for NGC 5466 are (m — M)y = 16.15 and 
E(B - V)=0.023 (Dotter et al. 2010). 

Four typical spectra are shown in Figure[2j to provide a 
sense for the quality and characteristics of the data. The 
top three spectra are a series in luminosity spanning the 
range of our data set, with the faintest at the top. The 
lowest spectrum is of star 1839, and shows the unusually 
strong CN and CH bands. 

2.1. The LRIS Comparison Sample 

Spectra were acquired for six RGB stars in NGC 5466 
by GHS in 2003, using the LRIS double spectrograph 
(Oke et al. 1995) at the W.M. Keck Observatory on 
Mauna Kea. All light was directed to the blue arm of the 
spectrograph by the use of a mirror in place of a dichroic. 
The 400/3400 grism was used, resulting in a typical res- 
olution of 7A and a pixel spacing of lA. Exposure times 
were between 420 and 600 s producing a mean S /N of 225 
pixel -1 in the wavelength range 4000A < A < 4100A. All 
reductions from flat field division through flux calibration 
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Fig. 3. — Radial velocities for all stars in our photometric sample, 
with a horizontal dashed line marking the mean velocity of 118.0 
kms" 1 . The horizontal solid line represents the mean cluster radial 
velocity values taken from Peterson & Latham (1986) and Pryor 
et al. (1991). The error on our mean velocity is 0.4 kms" 1 , and 
the rms is 17.0 kms -1 . This rms is 5% of a resolution element 
with VIRUS-P. The single inverted triangle represents star 1980, a 
strong-lined star that is also a radial velocity non-member. 



were carried out using the XIDL routines developed by 
J.X. Prochaska. 

3. RADIAL VELOCITY AND MEMBERSHIP 

For radial velocity determinations, data were cross- 
correlated with a template spectrum: for RGB and AGB 
stars, the template was the Hinkle Arcturus atlas (Hin- 
kle et al. 2000) convolved with a Gaussian to create 
an R = 850 spectrum. The HB star template was gen- 
erated from a synthetic spectrum with T ff= 6000 K, 
log(g)= 3.00 dex and [Fe/H]= -2.0 dex. 

Figure [3] shows the distribution of our velocity mea- 
surements. The weighted average is 118.0 ± 0.4 
kms _1 with an rms of 17.0 kms -1 . This rms is 5% of a 
resolution element with VIRUS-P. More customized tem- 
plate spectra may allow improvements in the rms scatter. 
The mean of literature values for the systemic velocity 
of NGC 5466 is 108 km s" 1 (Harris 1996), and solid hori- 
zontal lines in Figure 4 show velocities from Pryor et al. 
(1991) and Peterson & Latham (1986). The individual 
star radial velocities are listed in Table Q] In the fol- 
lowing analysis, spectra are individually shifted to rest 
wavelengths according to these radial velocities. 

4. BAND STRENGTHS AND BIMODALITY 

CN and CH band strengths are typically quantified 
with indices that measure the magnitude difference be- 
tween the integrated flux within the feature in question 
and the integrated flux in one or two nearby continuum 
regions, in the sense that more absorption in the feature 
produces a larger band-strength index. For CN band 
strength, we measure the index 5(3839), defined in Nor- 
ris et al. (1981) as 



5(3839) = -2.5 log 



f 3916 T 
J3883 IxCtA 



(1) 



. 2635 3176 727 1320 2891 



i S,(CH) - 0.029 +/- 0.005 



2635 3176 727 1320 2891 



Fig. 4. — Comparison of band indices measured from the Keck 
LRIS spectra discussed in Section 2.1 and those measured from 
VIRUS-P spectra, for the five stars in common to both data sets. 
Identification numbers are given at the bottom of each panel. 

The index we use for CH band strength is S2(CH), de- 
fined in Martell et al. (2008b) as 



S 2 (CH) = -2.5 log 



r 4317 t j\ 
J4297 -'agm 



f 4242 r r 4375 T 

J4212 + J4330 Ja " A 



(2) 



Band-strength indices tend to be designed to be effec- 
tive for a single stellar population, for example CH (4300) 
(Harbeck et al. 2003, designed for main-sequence stars 
in 47 Tucanae) and CH(G) (Lee 1999, designed for 
bright red giants in M3), and can be difficult to ap- 
ply outside the parameter space for which they were in- 
tended. S2(CH) was specifically designed to be sensitive 
to variations in carbon abundance in bright red giants 
(My < —1.0) across a wide range of metallicity, and 
therefore ought to be valid for the sample at hand. 

The errors on our band strengths have two compo- 
nents, errors related to the fluxing and errors related to 
the S/N. S/N errors were calculated by propagating the 
Poisson noise errors through the formulas for the band 
indices. To estimate the errors in our fluxing we tried 
several techniques to reflux the data including using syn- 
thetic spectra matched to each star's photometric tem- 
perature and using average flux distributions for color 
bins along the giant branch. The two sources of error 
were added in quadrature, and the total errors are re- 
ported in Tablets a s and a S2 for 5(3839) and S 2 (CH), 
respectively. For the March observing run the S/N in the 
52 (CH) band was quite high and the fluxing errors domi- 
nate while for the July observing run, where only a single 
dither set was taken for the final fourth field, the Poisson 
noise contribution plays a significant contribution to the 
total error. For 5(3839), the Poisson error tends to dom- 
inate the total error in all cases. Errors in radial velocity 
are typically 17 kms" 1 , or 5% of a pixel, as mentioned 
in Section 3, and do not have a significant effect on mea- 
sured band strengths. Shifting our wavelengths to the 
blue or red by 17 kms _1 changes the measured S 2 (CH) 
by 0.001 mag, which is small relative to the fluxing- and 
noise- related errors in S 2 (CH). 

Five of the stars in the LRIS data set are also in the 
VIRUS-P spectroscopic sample, and the CN and CH 



5 




□ ' 










□ □ ' 




□ 




□ 



SS(3839) 



Fig. 5. — CN band-strength index 5(3839) vs. absolute V mag- 
nitude, for all stars in our spectroscopic sample. The solid line is a 
least-squares fit, and we classify stars that lie above the line (shown 
as filled squares) as relatively CN-strong and stars that lie below 
the line (drawn as open squares) as relatively CN-weak. Small tri- 
angles denote stars with sufficient S/N for spectroscopic analysis 
but possible spectroscopic blending. The clear outlier, outlined 
with a large open star, is star 1839, noted in Figure 2 as a proba- 
ble CH star. The two points outlined with large open circles (stars 
1398 and 2483) lie to the blue of the RGB in Figure 1, and arc 
likely AGB stars. 



Fig. 6.— CN band-strength index 55(3839) vs. CH band- 
strength index S2{CH), for the spectroscopic sample. The symbols 
are the same as in Figure 5. The large crosses are centered at the 
mean values of (55(3839) and S 2 (CH) for the CN-weak and CN- 
strong groups, and their sizes are the standard deviations of the 
indices. CN and CH band strengths are typically anticorrelated 
in globular cluster stars, indicating anticorrelated carbon and ni- 
trogen abundances. The lack of strong anticorrelation here is the 
result of a smaller than typical range in CN band strength in NGC 
5466, and possibly a smaller range in [C/Fe] and [N/Fe] than in 
higher-metallicity globular clusters. 



band strengths measured from the LRIS spectra are quite 
similar to band strengths measured from the VIRUS-P 
data. Figure ffl shows A5(3839) and AS 2 (CH) for the 
stars in common between the two data sets, in the sense 
LRIS— VIRUS-P, and differences between the two data 
sets are quite small. The band strengths given in Table 
[2] come exclusively from the VIRUS-P data. 

Since the 3883A CN band strength can be used as a 
proxy for nitrogen abundance, and the 4320A CH band 
strength for carbon abundance, a comparison of CN and 
CH band strength in a globular cluster ought to re- 
veal anticorrelated abundance behavior. As can be seen 
in Figure El the CN band strength in our sample rises 
slightly with rising luminosity, but does not show a large 
range at fixed My. The clear outlier, with a large open 
star drawn around it, is star 1839, the star with the un- 
usually strong CN and CH features in Figure 2. It is 
excluded from the rest of our band strength and abun- 
dance analysis. 

The 5(3839) band strengths are correlated, to first or- 
der, with luminosity. To better distinguish stars that are 
CN-weak from stars that are CN-strong, we fit a linear re- 
lationship between 5(3839) and My and evaluated each 
star based on its vertical distance from this linear rela- 
tionship (55(3839)). We divide the sample into relatively 
CN-weak (negative 55(3839), open squares) and CN- 
strong (positive 55(3839), solid squares) stars. Unlike in 
studies of higher-metallicity globular clusters (see Norris 
et al. 1981 for a typical example), there is not a clear gap 
between the relatively CN-weak and CN-strong groups. 
In Figure [5] we compare (55(3839) and 5 2 (Cff) band 
strengths, with the same symbols (open/solid squares) 
representing CN class, and there is not a strong anticor- 
relation between CN class and CH band strength. The 
two large crosses are centered at the mean values for 
(55(3839) and S2{CH) for the two groups, with their size 
set to equal the standard deviations. They are not ob- 



viously offset from each other in S2(CH). This may be 
simply an effect of the extremely limited range in CN and 
CH band strengths permitted by the low overall metal- 
licity of NGC 5466. 

Another way to visualize the distribution of CN band 
strength in our data set is the generalized histogram, 
which is constructed by representing each star as a Gaus- 
sian in (55(3839) with a width equal to the measurement 
error erg, and then summing the individual Gaussians. 
The result is shown in Figure [JJ in which the solid curve 
is the generalized histogram for the full data set and 
the dashed curves are the generalized histograms cal- 
culated just for the relatively CN-weak and CN-strong 
groups. The upper panel of Figure [5] shows the single 
Gaussian that best fits the generalized histogram from 
Figure [7J and the lower panel, showing the residual be- 
tween the best-fit curve and the data, shows a clearly 
double-peaked shape. This double-peaked residual sug- 
gests the system is better fit by two Gaussians. 

Based on our analysis of the generalized histogram, we 
claim that NGC 5466 contains two closely spaced groups 
in CN band strength, a result of intrinsic ranges in car- 
bon and nitrogen abundance and low overall metallic- 
ity The lack of a gap in CN band strength between 
the two groups of stars in Figure [5] is seen here as a 
fairly small distance between the peaks of the CN-weak 
and CN-strong generalized histograms. Similarly metal- 
poor globular clusters such as M15 (Lee 2000) and M53 
(Martell et al. 2008a) also exhibit small separations be- 
tween their CN-weak and CN-strong groups and a small 
overall range in 5(3839) (among red giants in the three 
clusters (55(3839) has a standard deviation of 0.038 in 
NGC 5466, 0.056 in M53, and 0.037 in M15). This sug- 
gests that the range of carbon and nitrogen abundance 
in low-metallicity globular clusters may be smaller than 
is typical for more moderate-metallicity clusters. In ad- 
dition, since the measurement errors on 5(3839) control 
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Fig. 7. — Solid curve is a generalized histogram for the 
luminosity-independent CN band strength measure 55(3839) mea- 
sured for our full spectroscopic sample, and the dashed curves are 
generalized histograms for ^S^SSSE)) in the relatively CN-weak and 
CN-strong groups calculated independently. 

the widths of the individual Gaussians comprising the 
generalized histogram, they can have strong effects on 
our interpretation of it, with very large measurement er- 
rors potentially rendering two closely spaced populations 
indistinguishable from a single broad distribution. 

To understand the effects of measurement error on our 
generalized histogram, we performed a Monte Carlo test 
exploring how the size of our measurement errors af- 
fects the Kolmogorov-Smirnov coefficient p, which quan- 
tifies the probability that the CN-strong and CN-weak 
groups were drawn from the same initial population. We 
rescaled our true measurement errors ig by factors rang- 
ing from 0.5 to 10, then selected a new value for each data 
point from within a Gaussian distribution centered at the 
measured value of 5S*(3839), with a width of the rescaled 
measurement error. This reselection was done for the 
CN-normal and CN-strong groups identified in Figure [4] 
independently. The K-S p coefficient was then measured 
as an indicator of whether the two CN band strength 
groups containing reselected points with rescaled errors 
seemed likely to have been drawn from the same initial 
population. Reselections were done 10 4 times for each 
rescaling of the measurement error, and we find that p is 
strongly concentrated near zero when the rescaling fac- 
tor is relatively small: for a rescaling factor of 1, 99% 
of realizations have a p less than 0.1 and none have a p 
larger than 0.61, and 82% of realizations with a rescaling 
factor of 3 have a p less than 0.5. With a rescaling factor 
of 10, we find that p values are fairly evenly distributed 
between and 1, with 51% below p = 0.5. Based on 
this test, we are fairly confident that our two CN groups 
really are distinct, even if our errors in S*(3839) are un- 
derestimated by a factor of 3. 

5. EVOLUTION OF CARBON ABUNDANCE 

Carbon abundances were determined by spectral syn- 
thesis of the G-band region. For each star, a set of syn- 
thetic spectra were created with a single T c g, log(g), 
and [Fe/H], and a range of possible values for [C/Fe], 
smoothed to the resolution of the VIRUS-P data. The 
synthetic spectra were generated using a modified version 
of MOOG 2009 (Sneden 1973) which generates fluxed 
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Fig. 8. — Upper panel shows the generalized histogram from 
Figure 7 as a solid line, with a simple Gaussian fit to the distri- 
bution overplotted as a dotted line. The lower panel shows the 
difference between the two curves, and there is a clear two-peaked 
shape, which suggests that the data are best represented as two 
overlapping Gaussian distributions. 

spectra instead of normalized spectra. We used the 
PLEZ2000 grid of metal-poor model atmospheres (Plez 
2000). The temperatures and surface gravities were de- 
rived from the photometry in Table Q] using the prescrip- 
tions in Ramirez & Melendez (2005). We adopted [a/Fe] 
= 0.3, 12 C/ 13 C= 6 and [Fc/H] = -2.0 for all stars. The 
choice of the carbon isotope ratio does make a small 
difference to the derived total carbon abundance. At 
these metallicities, where the lines are reasonably weak, 
a change of 12 C/ 13 C from 6 to 50 would make about 0.06 
dex difference to [C/Fe] which is much smaller than the 
quoted errors. The synthetic spectrum yielding the resid- 
ual (in the sense observed - synthetic) that is closest to 
zero in the region of the CH G-band was chosen to have 
the correct carbon abundance. Table [2] lists ID numbers, 
VIRUS-P CN and CH band strengths, and carbon abun- 
dances along with their associated errors for the stars in 
the spectroscopic sample and for those that have good 
S/N but are possibly spectroscopically blended. 

Observational studies of CH band strengths and [C /Fe] 
abundances in red giants in globular clusters and in the 
field (e.g., Gratton et al. 2000; Smith & Martell 2003; 
Smith et al. 2005) have shown clearly that bright, low- 
mass red giants in all Galactic environments experience 
steady carbon depletion as they evolve along the giant 
branch. This depletion only occurs in stars brighter than 
the "bump" in the RGB luminosity function, and pro- 
ceeds at a rate that is metallicity-dependent. The homo- 
geneity and wide luminosity range of our spectroscopic 
sample makes it very useful for careful study of the rate 
and point of onset of deep mixing in NGC 5466. 

The first dredge-up process that occurs when a star 
moves on to the RGB (Iben 1965, 1968) is caused by the 
base of the convective envelope dropping in to small ra- 
dius, until hydrogen shell burning ignites, equilibrates, 
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Fig. 9. — Upper panel shows the CH band-strength index 
S2(CH) as a function of absolute V magnitude, and the usual 
globular cluster pattern of a plateau at low luminosity, followed by 
a decline in stars brighter than the "bump" in the RGB luminosity 
function, is apparent. The vertical line marks the position of the 
RGB bump at this metallicity. The lower panel shows our calcu- 
lated [C/Fe] abundances, which also exhibit a plateau followed by 
a steady evolutionary decline. As discussed in the text, the rate of 
carbon depletion measured in NGC 5466 is consistent with the low- 
metallicity carbon depletion rate found in Martell et al. (2008c). 



and restores the stability against convection of much of 
the stellar interior, driving the base of the convective 
zone to larger radius. This event homogenizes all ma- 
terial in the star outside the minimum radius reached 
by the convective envelope, and creates a discontinuity 
in the mean molecular weight \x at that radius. This 
fi discontinuity is destroyed when the hydrogen-burning 
shell, constantly moving outward in radius in search of 
fuel, encounters it. When that happens, the sudden 
change in the composition of the fuel available to the 
hydrogen-burning shell causes a brief loop in the star's 
path through the color-magnitude diagram: it becomes 
hotter and dimmer before recovering its equilibrium and 
continuing to evolve along the RGB. This evolutionary 
stutter causes a pile-up of stars in the luminosity func- 
tion of the RGB, known as the "RGB bump" (see, e.g., 
Fusi Pecci et al. 1990). The time at which a red giant's 
hydrogen-burning shell crosses the /x discontinuity is a 
unique function of stellar mass and metallicity, meaning 
that the absolute V magnitude of the RGB bump can be 
calculated empirically or theoretically for a given globu- 
lar cluster. 

This destruction of the \i barrier permits slow mix- 
ing between the surface and the hydrogen shell, which 
results in a continuous decline in surface carbon abun- 
dance and an associated decline in CH band strength 
starting at the RGB bump and continuing for as long as 
a star is on the giant branch (see, e.g., Palacios et al. 
2006; Denissenkov & VandenBerg 2003). The structure 
of the hydrogen-burning shell is also a function of stellar 
metallicity, and is more compressed in higher-metallicity 
stars (see, e.g., Sweigart & Mengel 1979 for a discus- 
sion of hydrogen-burning shell structure) . The more ra- 
dially extended structure of the hydrogen-burning shell 
in lower-metallicity stars allows deep mixing currents to 
carry material more efficiently from the stellar surface 
into the zone in which the CN cycle is operating, and re- 



turn carbon-depleted, nitrogen-enhanced material to the 
surface. The rate of carbon depletion at the surface, 
then, carries a fair amount of information about the in- 
ternal structure of the star, and also about the process 
that drives deep mixing currents (e.g., Eggleton et al. 
2008). 

The upper panel of Figure [H] shows CH band strength 
as a function of absolute V magnitude, and the expected 
pattern of a plateau followed by a steady decline with ris- 
ing luminosity can be seen. We calculate that the RGB 
bump ought to be at My = —0.41, using the relation 
My mp = 0.94[Fe/H]+1.5 found empirically in Fusi Pecci 
et al. 1990, and the vertical line in Figure|9]is at that lu- 
minosity. It is unclear whether our data completely agree 
with this prediction; given the sparseness of the data in 
the range —0.5 < My < 0.0, the corner between the 
plateau and the evolutionary decline could be anywhere 
within that range. 

The lower panel of Figure [9] shows our derived car- 
bon abundances as a function of absolute V magnitude, 
and [C/Fe] shows the same plateau, corner, and decline 
above the luminosity of the RGB bump (My ~ —0.4) as 
S^CH) does. The carbon abundance before the onset of 
deep mixing is roughly solar, and drops by a dex by the 
tip of the giant branch. As was demonstrated in Martell 
et al. (2008c), carbon depletion in globular cluster red gi- 
ants proceeds roughly twice as fast in the low-metallicity 
clusters M92 and M15 as it does in the higher-metallicity 
clusters M5 and NGC 6712. Using Yale-Yonsei evolution- 
ary tracks (Demarque et al. 2004), we determine the age 
difference between a star at My = —0.41 (the height of 
the RGB bump) and My = -1.5 to be 20 Myr, imply- 
ing a carbon depletion rate of roughly 34 dex Gyr -1 . 
This is in good agreement with the result of Martell 
et al. (2008c), which found a carbon depletion rate of 
20 — 30 dex Gyr -1 for globular clusters with metallicities 
between —2.0 and —2.5. 

6. SUMMARY 

We do not see strong star-to-star variations in CN or 
CH band strength in NGC 5466, but we do find that 
a sum of two Gaussians represents the generalized his- 
togram of CN band strength better than a single Gaus- 
sian curve. This suggests that the underlying 5(3839) 
distribution is likely to consist of two groups that are 
distinct but not widely separated. There is a very weak 
anticorrelation between CN and CH band indices, but it 
is muted even compared to similarly metal-poor clusters. 
It may be that NGC 5466 experienced weaker primordial 
enrichment than other Galactic globular clusters, and 
even other low-metallicity globular clusters. It would 
be interesting to investigate whether theoretical scenar- 
ios for primordial enrichment of globular clusters pre- 
dict varying enrichment efficiencies for particularly low- 
metallicity or low-mass clusters. 

Although the initial abundances in NGC 5466 may ex- 
hibit less of a range than abundances in typical globular 
clusters, the stellar evolution-driven abundance changes 
match well with what is expected from the literature. 
Surface carbon depletion, a result of slow circulation 
of material between the photosphere and the hydrogen- 
burning shell in bright red giants, begins near the 
"bump" in the RGB luminosity function and proceeds 
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at a rate consistent with the observations of Martell et 
al. (2008c) and the predictions of Denissenkov & Van- 
denBerG (2003). 

The data set presented here is particularly well suited 
to the study of deep mixing in the red giants in NGC 
5466, since the data extends to quite faint magnitudes. 
This is made possible by the multiplexing ability of 
VIRUS-P as an integral held spectrograph. We also show 
that VIRUS-P can be used to get velocities for stars to 
better than 17 kms _1 if good template matching is done. 
The future VIRUS instrument on the Hobby-Eberly Tele- 
scope will have comparable spectral resolution. However, 
instead of a single spectrograph it will have more than 
100 IFU spectrographs, and the fiber size on the sky will 
shrink from 4.1" to 1.5". This will make VIRUS a very 



capable instrument for observing larger, more crowded 
globular clusters. 
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t A ,r\rz.. A A c; q 
14:Uo:44.oo 


OO . Q 1 .1 c; 

2o:ol:lo.o 


15.56 


14.62 


13.55 


103.0 


5.6 


2629 


i/i ocr on 00 

14:U5:o9.2o 


0001 ion 

2o:ol:12.9 


14.89 


13.61 


12.27 


108.1 


7.0 


2642 


1 /i.nCi/io in 
14:Uo:4Z. iy 


00.01.10 1 
2o:ol: 12. 1 


16.76 


16.71 


16.63 


111.0 


13.0 


3176* 


1 a t~i tr on n/ ■ 

14:05:39. y6 


00 on 11 1 

2o:o0:ll.l 


15.61 


14.60 


13.50 


133.7 


13.0 


3306 


i/i ncr /ii no 

14:U5:41.Uo 


00 on a 

2o:29:4o.2 


15.90 


14.97 


13.92 


131.7 


6.0 


837 


1/t.nCiQn on 
14:Uo:oy.oU 


OO . o . c; O Q 

2o:oo:o2.o 


17.13 


16.36 


15.46 


121.2 


a 7 
b. / 


1398 


14:05:29.10 


28:32:47.3 


15.32 


14.48 


13.48 


123.2 


7.4 


1 Q/1 A 

lo44 


1 /i -nc;-Ofi 71 
14.Uo.2D. ( 1 


zo.oz.ol.D 


18.66 


18.20 


1 7 Q/1 

1 / .o4 


1 nQ 1 
lUo. 1 


/11V 

41. ( 


1106 


i/i ncr 00 tr 

14:05:32.58 


28:33:16.2 


18.89 


18.40 


17.63 


75.8 


17.1 


3387 


14:05:20.09 


28:29:24.1 


19.21 


18.66 


17.88 


132.8 


6.9 


2725 


14:05:22.37 


28:30:58.4 


18.82 


18.23 


17.41 


123.8 


29.8 


2957 


14:05:22.72 


28:30:33.9 


18.09 


17.64 


16.82 


94.9 


9.4 


3192 


14:05:22.77 


28:30:03.6 


19.46 


19.29 


18.78 


152.9 


22.5 


3078 


14:05:23.48 


28:30:21.4 


19.25 


18.67 


17.90 


144.6 


14.4 


946 


14:05:33.66 


28:33:36.3 


19.03 


18.46 


17.71 


108.8 


15.2 


968 


14:05:35.90 


28:33:33.7 


19.49 


18.88 


18.16 


113.8 


5.1 


888 


14:05:36.45 


28:33:45.2 


18.64 


18.08 


17.28 


115.3 


8.1 


3113 


14:05:26.33 


28:30:17.1 


18.96 
Sample 


18.35 
non-members 


17.58 


122.3 


16.3 


1980 


14:05:31.85 


28:32:02.6 


17.95 


17.21 


16.47 


68.9 


6.3 
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TABLE 2 

Results 



ID 


5(3839) 




S 2 (CH) 


o"52 


[C/Fc] 


CT [C/Fe] 








Main sample 








1839 


0.460 


0.027 


0.793 


0.038 


0.02 


0.18 


1398 


-0.185 


0.034 


0.603 


0.042 


-0.68 


0.18 


2483 


-0.199 


0.026 


0.609 


0.039 


-0.81 


0.14 


2629 


-0.093 


0.055 


0.510 


0.046 


-1.14 


0.26 


3176 


-0.028 


0.033 


0.622 


0.041 


-0.69 


0.13 


3306 


-0.160 


0.031 


0.663 


0.040 


-0.60 


0.43 


2635 


-0.144 


0.031 


0.596 


0.040 


-0.87 


0.17 


3127 


-0.040 


0.024 


0.522 


0.037 


-0.91 


0.13 


904 


-0.087 


0.023 


0.591 


0.038 


-0.80 


0.18 


727 


-0.135 


0.025 


0.648 


0.038 


-0.49 


0.14 


3071 


-0.166 


0.037 


0.587 


0.042 


-0.86 


0.13 


2891 


-0.110 


0.025 


0.636 


0.038 


-0.58 


0.11 


1320 


-0.139 


0.024 


0.630 


0.038 


-0.58 


0.11 


3305 


-0.144 


0.026 


0.658 


0.038 


-0.48 


0.12 


1497 


-0.102 


0.026 


0.645 


0.038 


-0.52 


0.11 


1534 


-0.166 


0.024 


0.648 


0.038 


-0.47 


0.11 


2776 


-0.145 


0.026 


0.664 


0.039 


-0.39 


0.11 


2821 


-0.209 


0.029 


0.688 


0.040 


-0.25 


0.12 




f» 1 1 R 
-U.11D 


n noo 
u.Uzy 


U.ooo 


U.U4U 


-U.oz 


n 11 
U. 11 


2544 


-0.218 


0.028 


0.656 


0.040 


-0.42 


0.13 


/Oi 


-0.206 


0.026 


0.709 


0.039 


-0.08 


0.12 


1/49 


-0.168 


0.039 


O.O/O 


0.045 


-0.11 


0.13 


826 


-0.138 


0.034 


0.685 


0.042 


-0.05 


0.10 


1946 


-0.185 


0.037 


0.701 


0.044 


-0.05 


0.11 


3401 


-0.203 


0.032 


0.685 


0.041 


-0.01 


0.14 


690 


-0.265 


0.054 


0.736 


0.054 


0.16 


0.15 


2753 


-0.239 


0.044 


0.709 


0.048 


-0.05 


0.21 


837 


-0.214 


0.031 


0.710 


0.041 


0.02 


0.13 


2957 


-0.238 


0.091 


0.660 

Sample with potential blends 


0.079 


0.12 


0.12 


2919 


-0.209 


0.042 


0.726 


0.046 


0.01 


0.12 


1351 


-0.188 


0.031 


0.712 


0.041 


-0.05 


0.14 


1075 


-0.220 


0.036 


0.704 


0.044 


0.05 


0.11 


1361 


-0.184 


0.053 


0.663 


0.053 


-0.03 


0.13 


1615 


-0.215 


0.036 


0.685 


0.043 


0.06 


0.10 


1024 


-0.191 


0.038 


0.659 


0.045 


-0.12 


0.15 


2481 


-0.174 


0.039 


0.560 


0.043 


-0.96 


0.11 


1972 


-0.088 


0.066 


0.631 


0.055 


-0.56 


0.11 



